Anomalous transport: a deterministic approach. 
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We introduce a cycle-expansion (fully deterministic) technique to compute the asymptotic behav- 
ior of arbitrary order transport moments. The theory is applied to different kinds of one-dimensional 
intermittent maps, and Lorentz gas with infinite horizon, confirming the typical appearance of phase 
transitions in the transport spectrum. 
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Generic dynamical systems are characterized by the 
coexistence of chaotic regions and regular structures, 
and typical trajectories present regular segments, due to 
sticking to the ordered component of the phase space, 
separated by erratic behavior, due to wanderings in the 
chaotic sea. Though ubiquitous, this mixed behavior still 
involves hard theoretical problems, as present theories 
are tailored to fit the two opposite paradigms of either 
integrable or fully chaotic systems. A particularly rele- 
vant feature associated to such weakly chaotic dynamics 
is anomalous transport: important features of the pro- 
cess are captured by the function u{q) that expresses the 
asymptotic growth of moments of arbitrary order: 



(1) 



ordinary diffusion yields u{q) = g/2, while anomalous 
transport often leads to a non trivial behavior, which 
cannot be encoded by a single exponent, but rather typi- 
cally exhibits a phase transition 1, 2]. Anomalous diffu- 
sion has been recognized as an ubiquitous phenomenon 
in recent years, from intermittency induced anomalous 
transport in one-dimensional maps [3] to the analysis of 
area-preserving maps in the presence of self-similar regu- 
lar structures around accelerator modes Q , from passive 
tracers dynamics in rotating flows 01, to charge carrier 
transport in amorphous semiconductors^ , to many oth- 
ers physically relevant contexts (see and references 
therein). In this paper we show how the transport ex- 
ponents u{q) may be computed in a crisply deterministic 
manner, by means of periodic orbit expansions '^,'^: we 
will then apply the formalism to various low dimensional 
settings, where analytic computations may be performed. 
These computations accomplish a twofold goal: besides 
checking the theory, they will also illustrate how subtle 
features of the underlying dynamical system are auto- 
matically included in the formalism. 
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For simplicity, we illustrate the analytic technique for 
the case of a one-dimensional map on the real line, even 
though the method is by no means limited to this con- 
text. The key property of the dynamical system under 
investigation is represented by the symmetry properties 



fi-x) = -fix) 



f{x + n) =n + fix) (2) 



for any n G N. These properties guarantee the absence 
of a net drift, as well as that the map on the real line is 
obtained by lifting a circle map f{0) (on the unit torus) 



/» = fm 



mod 1 



e T = [0, 1) 



(3) 



We may thus split the / evolution into a box integer plus 
a fractional part: Xn = Nn + On, where 



Nn+l = Nn 



■cr{f{0„)) 



,+1 = fiOn) (4) 



where (j{y) = [y] is the denotes the integer part. Trans- 
port properties are accounted for by the generating func- 
tion Gn{P), defined as 



GniP) ={e 



(5) 



(the average being taken among initial conditions), whose 
asymptotic behavior may be associated to the behavior 
of a generalized dynamical zeta function |3. ITol ITU . 
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(6) 



and ( is expressed as an infinite product over prime pe- 
riodic orbits of the torus map /: 



^(0)/3(^) 



n 



1 



lAJ 



(7) 



The quantities that enter the definition {TJ are the 
prime period Up of the orbit p, its instability Ap = 

nr=o ^ fiPi^p) integer factor a^, that accounts 

for the orbit's behavior once we unfold it on the real 
line. As a matter of fact, given any point Xp belonging 



2 




FIG. 1: A torus map with intermittent fixed points, and the 
corresponding hft. 



to p (that is f^p{xp) = Xp) we may either have that it 
is a periodic point of the hft / (i.e. f^^ixp) = Xp), or 
it might be a running mode, f^^ixp) = Xp -\- <jp, with 
cTp G Z (see If the zeta function Q has a simple 

zero z{P), then, from it fohows that the second mo- 
ment of the distribution grows hnearly in time (normal 
diffusion). This is generally the case with fully chaotic 
systems Isl [13. |l3l [13 : in the last few years it has been 
realized jld. lid. llTl Il8| that weakly chaotic systems (in 
particular one dimensional intermittent maps, or infinite 
horizon Lorentz gas models 19]) lead to more compli- 
cated analytic structure of the zeta function Q, which 
typically exhibit branch points. In view of the inverse 
Laplace formula © the modified analytic structure may 
induce anomalous behavior (nonlinear diffusion '^,'2^). 

We address here the problem of going beyond diffu- 
sion, and characterize the whole spectrum of transport 
exponents u{q) (^). From small /3 expansion of ©, we 
get 



^/c(^) = {{Xn-Xo)^)o 
k I r 

¥2^1 J, 



d(3 



gk 



ds 
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/3=0 
C(0)/3(^ ^) 



(8) 

/3=0 



The evaluation of the integral on the right hand side of 
(|HI) requires dealing with high order derivatives of a com- 
posite function: this is accomplished by making use of 
Faa di Bruno formula: 



d^H 



{L(t))-BAL(t)) (9) 



1 dL\ 



ki 



1 d^L 



(10) 



1\ dt J \n\ dt"^ 
^ = {/ci, . . . /cn} with ^ /ci = /c, ^i-/ci=n (11) 

When the analytic structure of the zeta function is 
known, from ^ we may single out the leading singular- 
ity in the logarithmic derivative, and then estimate the 
asymptotic behavior of 0, for instance b y em ploying; 
Tauberian theorems for Laplace transforms 
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We will apply the technique to two classes of one di- 
mensional maps, where deviations from fully chaotic be- 
havior are provided by marginal fixed points, of inter- 
mittent type 22]. The torus map is shown in fig. QJ: 
it consists of three branches, the central one being hy- 
perbolic (with constant slope), while the other two in- 
clude a marginal fixed point: to unfold it on the real line 
it is sufficient to assign jumping numbers to branches. 
For instance we get the lift of fig. (^) once we assign 
cr_ = —1 to the left branch, (Jc = to the central one 
and (j+ = +1 to the right branch. The map is actually 
taken as a straightforward generalization of the Gaspard- 
Wang piecewise linear approximation Q of Pomeau- 
Manneville map: the intermittent behavior is determined 
by dynamics near the pair of parabolic fixed points, which 
is accounted for by the intermittency exponent 7 > 1, 



(the map goes like x^+i 



near the origin, with 



an analogous behavior at the twin fixed point in x = 1). 
The zeta function can then be written as 



^ k 

k=i 



where a and b are fixed by specifying the central region 
slope and the normalization condition C(o)o(-'-) ~ ^' 



a = 1/(7-1). 



(13) 



The appearance of the Bose function 5'^(^) = Xl^^i JF 
is due to sequences of orbits coming closer and closer to 
the marginal fixed points: their stability increases only 
polynomially with the period 13, 23], a clear signature of 
local deviation from typical hyperbolic behavior (which 
is ruled by exponential instability growth). To estimate 
the various contributions in Q we remind the behavior 
as 2; 1~ 



(1 - z)^-^ 
ln(l - z) 

C(2)+C2(l- 



z)\n{l 



D^{l-z) 



/i < 1 
/i = l 

/iG(l,2)(14) 
11 = 2 
/i > 2 



and moreover take into account that 



=0 L 







i odd 



zga+i-i{z) i even 



(15) 



Now take a generic term in Q: and denote it by P/ci.../c^ 
we have in view of (|15|) 



ki ...kn 



{z)Y 



V7 



(16) 

where the picks up the contributions from the prod- 
uct of Bose functions, and all j must be even, due to (|15|) . 
First we consider the case a G (0, 1), which corresponds 
to 7 > 2: we have V^^ ~ (1 — z)^", that, together 



3 



with , implies that the dominant singularity is of the 
form 



where p is determined by 
p — sup (ka - 



(1-z)^ 



-a)kj) =n 



(17) 



(18) 



Once plugged into 
which means that 
ballistic behavior 
method applies). 
a > 1: since the 
zero we get 
in P+ modify the 
sufficiently high j 



|H|) this leads to the estimate u{q) = 
the whole set of moments is ruled by 
(at least for even exponents, where the 
We now turn to the more subtle case 
dynamical zeta function has a simple 
~ (1 — z)^, while the terms appearing 
singular behavior near z = 1 only for 



ga+i-j{z) 



(1 - z)"-^' 

C(a + l-j) 



J > a 
j <a 



(19) 



If all {j} are less than a then the singularity is deter- 
mined by V~: keeping in mind that the highest k value 
is achieved by choosing j = 2 and /c2 = n/2, we get, by 
proceeding as before 



q < a 



(20) 



When q exceeds a we have to take into account possible 
additional singularities in and thus we get 



where p is determined by 



{i-zY 



(21) 



p= sup {k+Y,{j-a)k,) = \ 1^1^ _^ n>2(a-l) 



L -L nj J>Ot 

which, once we take (PU)) into account, yields 



n < 2(a - 1) 
(22) 



r ^/2 



q<2{a 
a q > 2{a 



1) 
1) 



(23) 



The set of exponents thus has a nontrivial structure, 
characterized by a sort of phase transition for = 2(a — 
1), a rather universal feature of many systems exhibiting 
anomalous transport 0. We notice that the parameter 
ruling the presence of a phase transition (and the explicit 
form of the spectrum) is a, the exponent describing the 
polynomial instability growth of periodic orbits coming 
closer and closer to the marginal fixed point, and thus 
describing the sticking to the regular part of the phase 
space: in the present example a = (7 — 1)~^, and thus 
the sticking exponent is easily connected to the intermit- 
tency index. 

We now discuss a further 1-d example: first we re- 
call that the former torus map has non trivial ergodic 



f(x)0 




FIG. 2: The intermittent map with constant invariant mea- 
sure. 



' q ' 



FIG. 3: Spectrum of the transport moments for the map (f24|l 
with 7 = 1.5: the best fit on numerical data is ^ = 0.50x-h0.04 
for the dotted line, and y = 0.98x + 1.82 for the full line. 



properties: an absolutely continuous invariant measure 
only exists for a > 1 (see for instance 24J and refer- 
ences therein): the ergodic behavior is much more com- 
plex when a < 1. In the new example 1] while the 
functional form of the map near marginal fixed points is 
identical to the former case, ergodic properties are com- 
pletely different (the measure is not singular). The torus 
map, again dependent on an inter mittency parameter 7, 
is implicitly defined on T = [—1,1) in the following way 



0<x < 1/(27) 
,7 (24) 
fix) 1/(27) < ^ < 1 



for negative values of x the map is defined as f{—x) = 
—fix) (cfr (|21)): see fig. Q. The peculiar ergodic fea- 
tures of the map, namely the existence of a constant in- 
variant measure for any value of 7 arise from the property 
Xly=/-i(a;) ^/f{y) = 1 Also in this case we can easily 
identify families of orbits coming closer and closer to the 
marginal fixed points, but evaluating their instability re- 
quires some care, as the slope in the chaotic region is 
not bounded from above. A piecewise linear approxi- 
mation (in the same spirit as [2^) is still possible [25j . 
but we must put particular attention on matching the 
summation property : the corresponding dynamical zeta 
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function is 

oo ^ 

C(o)/j(^) = 1 - C(« + 2)^ ^ ^ cosh(/3fc) (25) 

where again a = 1/(7 — 1). Thus the continuity of the 
measure deeply modifies the relationship between the in- 
termittency exponent and the instabilities of periodic or- 
bits shadowing the marginal fixed points: by repeating 
the steps of our former calculation we get that we always 
get a phase transition with 

\^ q — a q > 2a ^ ^ 

which may also be checked numerically (see fig. (|2|). 

Our last example shows how the proposed technique 
may be applied to higher dimensional systems, where a 
detailed layout of the full symbolic dynamics is beyond 
our present understanding. We consider the Lorentz gas 
with infinite horizon (square lattice of circular scatterers 
where a test particle freely moves, colliding elastically 
with the disks). The mechanism that in this case leads 
to deviation from normal behavior is due to the possibil- 
ity of arbitrarily long flights of free (collisionless) motion 
along corridors. In this example the reduced dynamics 
is that of a Sinai billiard: while a complete description 



of the full set of periodic orbits of such a system is not 
known yet a a family of periodic orbits approaching closer 
and closer the infinite free fiights is singled out: it consists 
(for the unfolded billiard) to orbits that jump n lattice 
spacings between collisions: their number and instability 
can be evaluated by geometric arguments 19j, yielding an 
istability that grows polynomially with n with a power 
3: so their role is analogous to modes leading to Bose 
functions in (|T2j) with a = 2. In view of (|23j) this leads 
to u{q) = q/2 for q < 2 and iy{q) = g — 1 for g > 2, like 
recently suggested in [2^: a careful analysis, based on 
([T^ moreover gives (72{n) ~ nlnn (see 27J), so that our 
method is also capable of picking up logarithmic correc- 
tions to the dominating power-law behavior. 

We have proposed a crisply deterministic technique to 
investigate the full spectrum of transport exponent for 
chaotic systems: in particular this method is capable 
of explaining the different phase transitions that possi- 
bly arise, without any direct information on the invari- 
ant measure: the only information to be plugged is local 
growth of instabilities near the marginal structures, to- 
gether with the appropriate jumping factor. 
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